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Liquid crystals in two dimensions undergo a first-order isotropic-to-quasi-nematic transition, pro- 
vided the particle interactions are sufficiently "sharp and narrow". This implies phase coexistence 
between isotropic and quasi-nematic domains, separated by interfaces. The corresponding line ten- 
sion is determined, and shown to be very small, giving rise to strong interface fluctuations. When 
the interactions are no longer "sharp and narrow", the transition becomes continuous, with non- 
universal critical behavior obeying hyperscaling, and approximately resembling the two-dimensional 
Potts model. 
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Fluids consisting of elongated molecules (rods, nee- 
dles) can form nematic phases. In the nematic phase, 
the molecules are aligned, which distinguishes it from 
the isotropic phase, where the molecular orientations are 
random. The nematic phase is the standard example of 
a liquid crystal. Consequently, its properties are of fun- 
damental importance. The nematic phase can be pre- 
pared in several ways. In thermotropic liquid crystals, 
the isotropic-to-nematic (IN) transition is temperature 
driven: starting in the high temperature isotropic phase, 
the nematic phase is reached by lowering the temper- 
ature. The Lebwohl-Lasher (LL) model l| provides a 
convenient theoretical framework. In this model, a unit 
vector di is assigned to each site i of a lattice, where 
di represents the orientation of the molecule at the i-th 
lattice site. The molecules interact with Hamiltonian 

H = -eJ2\di-d^, (1) 

where the sum is over nearest neighbors (factors of k^T 
are absorbed in the coupling constant e, with T the tem- 
perature, and &b the Boltzmann constant). Note that 
Eq.([T]) is invariant under inversion, di — > — di, which is 
characteristic of liquid crystals. In the original paper, 
Eq. ([1]) is studied in d — 3 dimensions with p = 2 [1[ . In 
this case, at low enough temperature, Eq.([T]) undergoes a 
first-order transition from an isotropic phase to a nematic 
phase. Nematic phases also occur in lyotropic systems, 
where density drives the IN transition. In a seminal pa- 
per 0, Onsager showed that infinitely slender rods in 
d — 3 dimensions also undergo a first-order transition 
from an isotropic to a nematic phase, at sufficiently high 
density. 

The IN transition in d = 3 dimensions is well under- 
stood. In contrast, the two-dimensional case, which will 
be the topic of the present investigation, is more con- 
troversial. More precisely, we consider Eq.Q in d = 2 
spatial dimensions, using two-component unit vectors dj. 
For p = 2, Eq.JT]) then becomes the XY model @, and 
a nematic phase with true long-range order is ruled out 



by the Mermin- Wagner (MW) theorem 0]. The two- 
dimensional XY model and its variants were thought 
to be without a phase transition for a long time, un- 
til Kosterlitz and Thouless (KT) proved that a phase 
transition does occur, and clarified its topological nature 
[5(. The KT transition is one from a (high tempera- 
ture) isotropic phase, with exponential decay of the an- 
gular correlations, to a (low temperature) quasi-nematic 
phase, with power-law decay of the correlations. In the 
XY model, the KT transition is continuous. By lower- 
ing the temperature, starting in the isotropic phase, the 
correlation length grows exponentially, until it diverges 
at the transition temperature, where the quasi-nematic 
phase sets in. Since the quasi-nematic phase has infinite 
correlation length, it is a critical phase. Consequently, 
the order parameter A and the susceptibility x scale as 

A oc L-M", X oc L 7/l/ , (2) 

with L the system size; (3, 7, and v are the critical ex- 
ponents of the order parameter, susceptibility, and cor- 
relation length, respectively. Since the correlation length 
diverges exponentially, the exponents themselves are un- 
defined. However, exponent ratios are still defined [1, @], 
via Eq.©. 

Since the XY model and the LL model are similar, the 
IN transition in two-dimensional liquid crystals is often 
assumed to be of the conventional KT type. Clearly, for 
Eq.(JT]) with p = 2 this is justified. If one then accepts 
a well-defined universality class for two-dimensional liq- 
uid crystals, it seems reasonable that Eq.|T|) for arbitrary 
p > 2, and indeed also lyotropic liquid crystals, are all 
qualitatively similar. The purpose of this Letter is to 
demonstrate that the IN transition in two dimensions is 
far more subtle. Our results are inspired by generalized 
XY models in d = 2 [7J, for which it has been proved 
that the KT transition can become first order Q. Since 
Eq.([T]) can be mapped onto the generalized XY model, 
using (cosx) 2p — 2~ v (1 + cos(2a;)) p , the consequences 
of this result should be relevant for liquid crystals as 
well. Indeed, for liquid crystals interacting via Eq.JI]) 
with large p, a first-order transition is found, including 
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FIG. 1: Evidence of a first-order phase transition in the ther- 
motropic liquid crystal of Eq.JTJ with p = 1000 at e ~ 2.5. (a) 
Scaled and shifted logarithm of P{E) for L — 10 (solid curve) 
and L = 15 (dashed curve); A reflects the line tension of the 
IN interface, (b) Typical snapshot obtained in the coexistence 
region. 



a phase coexistence region, characterized by a finite line 
tension. For smaller p, non-universal critical behavior is 
found, with exponent ratios that vary continuously with 
p, while obeying hyperscaling. Interestingly, the varia- 
tion of the exponent ratios we observe is qualitatively 
similar to that of the d — 2 Potts model |9j. Finally, 
we consider an off-lattice liquid crystal, for which similar 
non-universal critical behavior is found. 

We begin our investigation by considering Eq.([T|) in 
the limit of large p. In this case, the nearest neigh- 
bor interaction is extremely "sharp and narrow" . In 
other words, two neighboring molecules lower the energy 
only when they are closely aligned; otherwise, the inter- 
action quickly vanishes, which mimics the Kronecker-(5 
Hamiltonian of the Potts model [9(. In fact, for large p, 
Eq.([T]) qualitatively resembles the (/-state Potts model, 
with q oc p 1 / 2 [7]. The g-state Potts model in d = 2 
exhibits a first-order phase transition when q > 4, and 
so one may hope to see a similar transition in Eq.|T|) 
when p becomes large. To this end, we have Monte Carlo 
simulated Eq.([T]) with p = 1000, on a periodic L x L 
lattice. The simulations are performed using standard 
single-particle Metropolis moves, combined with a biased 
sampling scheme [T(|, and histogram reweighting |TT1 |. 
Evidence of a first-order transition is obtained from the 
distribution P{E), defined as the probability to observe 
the energy E during the simulation [T3| , For e ~ 2.5 
and p = 1000, we find that P(E) becomes bimodal, see 
Fig. [TJa) . The peak at low energy (I) reflects the quasi- 
nematic phase, the peak at high energy (II) the isotropic 
phase, and the region in between the peaks corresponds 
to phase coexistence (the latter is characteristic of first- 
order transitions). Simulation snapshots obtained in the 
region where P{E) attains its minimum strikingly con- 



firm phase coexistence, see Fig.[ljb). Here, a rectangular 
lattice was used, such that the interfaces form parallel to 
the short edge of the lattice, since this minimizes the 
total amount of interface in the system. Note that the 
coexistence is between an isotropic and a quasi-nematic 
phase: both phases lack long-range order in the thermo- 
dynamic limit. The decay of nematic order with system 
size in the quasi-nematic phase may, however, be very 
slow, giving the impression of a true nematic phase [13| • 
The interfaces in Fig. QJb) are not flat, and appear to be 
decorated with capillary waves. Our results even allow 
for an estimate of the line tension A between the coexist- 
ing domains. To this end, note that we have plotted the 
logarithm of P{E) divided by 2L in Fig. QJa), with the 
minimum between the peaks shifted to zero. The height 
of the peaks then reflects the line tension [3] . The re- 
sults of both system sizes are remarkably consistent, and 
yield A ~ 0.3 k^,T per lattice spacing. 

Next, we investigate what happens when the interac- 
tion of Eq. |T]) is no longer "sharp and narrow" . Clearly, 
for p — 2 in Eq.{T]), critical behavior of the XY model 
should be detected. Therefore, somewhere in the inter- 
val 2 < p < 1000, a crossover to first-order behavior must 
take place. We have performed additional simulations of 
Eq.([T]), and find that for p < 50, bimodal energy distri- 
butions P(E) can no longer be identified. In fact, the 
case p — 50 is near the borderline: for small systems, 
bimodal energy distributions do occur, at e w 1.86, but 
the free energy barrier AF = 2LX decreases with L, and 
so the bimodal structure does not survive the thermo- 
dynamic limit. Following Ref. ITU we conclude that the 
transition in Eq.{T]) is no longer first-order when p < 50. 
Interestingly, a similar phenomenon also occurs in the 
two-dimensional g-state Potts model [9]. Here, q = 4 
is the borderline case: when q < 4, the Potts model no 
longer exhibits a first-order phase transition. Instead, the 
transition becomes critical, with critical exponents that 
depend on q. If we accept that p in Eq.lJTJ) is analogous 
to the number of Potts states q, then p = 50 roughly 
corresponds to q = 4. It then becomes of interest to 
investigate the critical behavior of Eq.|T|) in the regime 
p < 50. If the analogy to the Potts model remains valid, 
critical exponents that depend on p are to be expected. 

To this end, we now consider Eq.JT]) using p = 10. 
Since P(E) is no longer bimodal, a different quantity 
must be used to locate the phase transition. For liquid 
crystals, a natural choice is the nematic order parameter 
S 1 , defined as the maximum eigenvalue of the orienta- 
tional tensor Q a p — YliLi {^diadip — 5 a p), with di a the 
a component (a = x, y) of the orientation di of molecule 
i, 5 a f3 the Kronecker delta, and N the total number of 
particles. For disordered phases, S will be small; for or- 
dered phases with strong alignment, S will be larger. We 
simulate Eq.((TJ) as before, and measure the distribution 
P(S), defined as the probability to observe the nematic 
order parameter S. First evidence of a phase transition is 
provided in Fig. (Ha). Shown is W = \nP(S) for several 
values of e. The striking feature of Fig. HJa) is the for- 
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FIG. 2: Nematic order parameter distributions W = In P(S) 
of Eq.{T]) for p = L = 10. (a) W in the absence of JS, for 
several e, including the "common tangent" construction for 
e = 2.5. (b) Bimodal form of W for e = 2.5, in the presence 
of JS, including the definition of AK and A. 
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FIG. 3: Finite size scaling analysis of Eq.{TJ with p = 10. (a) 
Ui as a function of e for several system sizes L. The intersec- 
tion point yields an estimate of e cr . (b) AK as a function of 
1/L, for e = 1.70; 1.66; 1.64; 1.62; 1.60 (top to bottom). 



mation of a "kink" in W, when e is sufficiently large. It 
then becomes possible to perform a "common tangent" 
construction. We thus find a rather special phase transi- 
tion, characterized by a change in the shape of W. The 
transition is one from a high-temperature (low e) phase, 
where "common tangents" in W do not occur, to a low- 
temperature (high e) phase, where they do. We note in 
passing that the probability distribution of the magne- 
tization in the two-dimensional XY model shows similar 
behavior 0]. In other words, the formation of a "kink" 
does not imply long-range order in Eq.|T]), which indeed 
is ruled out by the MW theorem [4(. 

To accurately locate the value of e above which the 
"kink" in W begins to form, it is convenient to add a 
term JS to the Hamiltonian, such that Eq.((T|) becomes 
H = M« ' dj\ p + JS, with S the nematic order 

parameter, and coupling constant J. The effect of J > 



FIG. 4: Critical properties of Eq.{TJ as a function of p. 



is to penalize the nematic phase. Provided a "kink" is 
present, W can be cast into bimodal form. Fig. [2fb) 
gives an example for e = 2.5. Here, J was tuned us- 
ing Pj^q(S) = Pj = q(S) exp(— JS), such that the two 
peaks in P(S) were of equal area The magnitude 

of the "kink" now shows up as a barrier, marked AK in 
Fig-IHb). To locate the transition, we use the finite size 
scaling approach of Ref. [TH To this end, AK is mea- 
sured as a function of e and L. The result is shown in 
Fig. EJb). For small e, AK decreases with system size, 
implying that the "kink" does not survive the thermo- 
dynamic limit. For large e, AK increases with system 
size, in which case the "kink" survives. For interme- 
diate e, AK remains roughly constant, implying that 
the transition from the low-e "kink-less" phase, to the 
high-e "kink" phase, passes through a critical point [l5l ]. 
at e ss 1.64. Further confirmation is obtained from the 
Binder cumulant U4 = (m 2 ) / (m 4 ) of the bimodal dis- 
tributions, with m = S — (S) , see Fig. [SJa) . At a critical 
point, the cumulant becomes L-independent [13] ■ The 
sharp intersection point in Fig. (3ja) provides strong ev- 
idence that Eq.([T]) indeed becomes critical. For p = 10, 
criticality is obtained at e cr ~ 1.637, in excellent agree- 
ment with the previous estimate. 

Having established the critical value of e, critical ex- 
ponent ratios can be measured in several ways. For ex- 
ample, j3/v is obtained by fitting the decay of the order 
parameter A at e cr to Eq.([2]), with A defined in Fig.[2jb). 
Similarly, 7/^ follows from the finite size behavior of the 
susceptibility X = ((S 2 ) - (S) 2 )/L 2 at e cr , see Eq.@. 
In addition, exponent ratios can be obtained using the 
method of Loison [18| . We find that both techniques are 
remarkably consistent: for p = 10 in Eq.([T]), we obtain 
j3/v w 0.175 and j/u w 1.645. Note that these ratios 
strikingly obey the hyperscaling relation r y/is + 2(3/i' = d. 
By repeating the above analysis for different values of p, 
the result of Fig. |4] is obtained. Shown are e cr (a), 0/u 
(b), and 7/1/ (c), as a function of p. We first note that the 
exponent ratios for p ^ 10 also obey hyperscaling. This 
result is important because it demonstrates the consis- 



4 



tency of our data, and provides additional confirmation 
that Eq.(TT]), for small p, indeed becomes critical at the 
transition point. An even more striking feature is that 
the exponent ratios depend on p. Such non-universal 
critical behavior may seem surprising, but has been ob- 
served before in different systems [19|, [20] • Indeed, based 
on the analogy to the Potts model, p-dependent criti- 
cal behavior was already anticipated. In fact, the non- 
monotonic variation of the exponent ratios we observe in 
Fig. 01 is also characteristic of the g-state Potts model. 
In this case, the exponent ratios assume their extrema 
at q k 3.33. Our results thus suggest that Eq.([T]) for 
p ~ 10 — 20, roughly corresponds to a q ~ 3 Potts model. 
The actual exponent values of the Potts model, however, 
do not agree with those of Fig. |U In other words, the 
analogy between Eq. (JTJ) and the Potts model is not ex- 
act. Note that the rotation symmetry of Eq.JT]) is not 
the permutation symmetry of the Potts model. An ex- 
act correspondence is therefore not to be expected. This 
is also manifested by the behavior of e cr , see Fig. Ufa). 
Whereas the increase of e cr with p for large p is consistent 
with the Potts model, the decrease at small p is not. In 
fact, for p — 2 in Eq.JT]), we expect XY critical behavior 
to occur. As Fig. [4jc]demonstrates, 7/1/ is indeed close 
to the XY value 7/4 |3|] in that case. 

We now consider an off-lattice liquid crystal, namely a 
fluid of soft rods. The rods are defined as rectangles, of 
length I and width w, capped at each end by a semi-circle 
of diameter w (we set l/w = 16, and I will be the unit of 
length) . The rods interact via a repulsive pair potential, 
whereby rod overlap is penalized with energy cost k = 2. 
The rods are simulated in the grand-canonical ensem- 
ble, i.e. at constant temperature T, chemical potential 
/Lt, and system area A, while the number of particles N 
fluctuates. We use a simulation square of size L x L, 
with periodic boundary conditions. While for Eq.|T]) a 
phase transition occurs above a certain e, here that role 
is played by /1. This also implies that the analogue of en- 
ergy in Eq. (TT]), becomes the number of particles N, since 
\x couples to N in the grand-canonical ensemble. The 
probability distributions P(S) and P{N) were measured 
for various fi. We find that bimodal distributions P(N) 
do not occur, strongly suggesting that a first-order tran- 



sition is absent. In other words, soft rods resemble Eq.|T|) 
in the limit of small p. Indeed, we find that In P(S) de- 
velops a "kink" when /x exceeds a critical value fi CT . For 
/i > /i cr , bimodal distributions lnP(S) can be realized, as 
in Fig.[2fb), by tuning J. Finite size scaling confirms that 
the system becomes critical at the point where the "kink" 
first appears. For soft rods, /x cr w 1.985, fl/v s=s 0.17, and 
j/v w 1.65 are obtained, consistent with hyperscaling. 
By comparing to Fig. [4] we conclude that soft rods re- 
semble Eq.© with p ~ 10 - 20. 

In summary, for a lattice liquid crystal with "sharp 
and narrow" interactions, the existence of a first-order 
transition was shown, including an estimate of the line 
tension between the coexisting domains. The line ten- 
sion is small, giving rise to strong interface fluctuations. 
When the interaction is no longer "sharp and narrow" , 
the transition becomes continuous. Liquid crystals then 
show non-universal critical behavior, with exponent ra- 
tios that depend on the "sharpness" of the interaction, 
but that do obey hyperscaling. In addition, the behav- 
ior of the exponent ratios follows the Potts trend. For 
lattice liquid crystals, the transition type can be selected 
using the parameter p in Eq.([T]). For off-lattice systems, 
such a parameter is not so easily identified. In soft rods, 
the interaction is clearly not "sharp and narrow" enough 
to induce a first-order transition. It remains of inter- 
est to identify off-lattice interactions that do facilitate 
first-order transitions in two-dimensional liquid crystals. 
Possibly, this could be achieved using Gay-Berne type 
potentials [2l|. A different application could be to use 
the present simulation methodology to study melting in 
two dimensions. Such investigations are a topic for future 
work. 
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